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ABSTRACT 

We report the detection of the Paa emission line in the z = 2.515 galaxy SMM J 163554.2+66 1225 using 
Spitzer spectroscopy. SMM J 163554.2+66 1225 is a sub-millimeter-selected infrared (IR)-luminous galaxy 
maintaining a high star-formation rate (SFR), with no evidence of an AGN from optical or infrared spec- 
troscopy, nor X-ray emission. This galaxy is lensed gravitationally by the cluster Abell 2218, making it ac- 
cessible to Spitzer spectroscopy. We measure a line luminosity, L(Paa) = (2.05 ±0.33) x 10 42 erg s" 1 , cor- 
rected for gravitational lensing. Comparing the Ha and Paa luminosities, we derive a nebular extinction, 
A(V) = 3.6 ± 0.4 mag. The dust-corrected luminosity, L(Paa) = (2.57 ± 0.43) x 10 42 erg s" 1 , corresponds to 
an ionization rate, £>o = (1-6 ±0.3) x 10 55 7 s -1 . The instantaneous SFR is ip = 171 ±28 M Q yr -1 , assum- 
ing a Salpeter-like initial mass function from 0.1 to 100 Mq yr -1 . The total IR luminosity derived using 
70, 450, and 850 pm data is Lir = (5- 10) x 10 11 Lq, corrected for gravitational lensing. This corresponds to 
ip = 90- 180 M© yr" 1 , where the upper range is consistent with that derived from the Paa luminosity. While the 
L(8/im)/L(Paa) ratio is consistent with the extrapolated relation observed in local galaxies and star-forming 
regions, the rest-frame 24 pm luminosity is significantly lower with respect to local galaxies of comparable 
Paa luminosity. Thus, SMM J163554. 2+661225 arguably lacks a warmer dust component (T D ~ 70 K), which 
is associated with deeply embedded star formation, and which contrasts with local galaxies with comparable 
SFRs. Rather, the starburst in SMM J163554. 2+661225 is consistent with star-forming local galaxies with 
intrinsic luminosities, Lir w 10 10 Lq, but "scaled-up" by a factor of ~10-100. 

Subject headings: infrared: galaxies — galaxies: formation — galaxies: high-redshift — galaxies:individual 
(SMM J163554.2+661225) — galaxies: starburst 



1. INTRODUCTION 

There has been a growing industry of deep, multiwave- 
length surveys, devoted to studying star formation and evo- 
lution of galaxies. Studies of data from these surveys have 
concluded that the global star-formation rate (S FR) den- 
sity r eached a maximum between z~ 1.5-3 (e.g., iHopkinsI 
12004 and references therein), durin g the period of rapid 
growth in the stellar mass density (e.g. iDickinson et"ai]|2003t 

1 This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion laboratory, California In- 
stitute of Technology, under NASA contract 1407. 

- papovich@physics.tamu.edu 
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iRudnicket aTJ 120031 12006b . However, the study of SFRs in 
these multiwavelength datasets requires estimates of the SFR 
using available observables. Common SFR indicators for 
high-redshift galaxies are the rest-frame UV luminosity, es- 
pecia lly at z > 1 when t he UV shifts into the optical bands 
(e.g. iMadau et al.lll996t iSteidel et alj|1999t iGiavalisco et al l 



20041: iBouwens et al.l (20061: iSawicki & Thompson! I200i. 
Reddv et al.H2.Q08l) , the far-IR lumino sity inferred from sub 



mm emission (e.g.. iBlain et alj|2002l). the mid-IR luminos 
ity from Spitzer/24 /xm (e.g., Perez-Gonzalez et"aTJ 1 2005 
Caputi et al l 120061: iPapovich et alj 120061: iReddv etafl 2006 



Webb et al.H 2006t iDaddi et alJ|2007t iFranx et all 120081). and 



the Ha-emission l ine (e .g., Erb et al.l 120061: iBouche et al.1 
l2007t iKriek et al.l 120081) . Each of these SFR indi- 
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cators has inherent unc ertainties (e.g., iKennicutH 119981 : 
iHook ins & Beacom 2006, and references therein), and some 
that are unique to observations at these redshifts (e.g., un- 
certainties in the luminosity-temperature dependence on the 
sub-mm-LiR conversion, Chapman et al. 2005; Pope et al. 
2006; template incompleteness and photometric redshift un- 
certainties on the 24 /mi-LiR conversion, Papovich et al. 2006, 
Reddy et al. 2006). Nevertheless, several studies have com- 
pared SFRs derived from the UV, 24 fim, far-IR, and sub- 
mm (and less-understood SFR tracers, such as radio and X- 
ray emission), concluding broadly that whi le they are consis- 
tent st atistically, the re is la rge scatter (e.g..|Daddi et alJl2005l 
2007| iReddv et al.l [2006t iPope et al.1 120061: iPapovich et alj 
20071) . There have been few comparisons of direct, robust 



SFR indicators in individual galaxies at high redshifts (e.g., 
Siana et al. 2008; references above). 

SFR indicators ar e normally calibr ated against H I recom- 
bination lines (e.g.. lKennicuttl[T998l) . The primary advan- 
tage of H I recombination lines is that they effectively re- 
emit the stellar luminosity from hydrogen-ionizing photons. 
They are strong in star-forming regions, and they trace di- 
rectly emission by massive stars. The physics of H I emis- 
sion in star-forming regions is relatively well understood (Os- 
terbrock 1989), and they exhibit only a weak dependence 
on electron density and temperature. Of the H I lines, Paa 
at A = 1.8751 /im is an ideal practical SFR indicator. This 
line is exceptionally strong in star-forming regions (1 M Q 
yr" 1 corresponds to L[Paa] = 1.48 x 10 40 erg s" 1 ; Alonso- 
Herrero et al. 2006). Furthermore, while recombination lines 
are subject to strong attenuation, Paa suffers minimal extinc- 
tion due to its longer wavelength. Recent studies of star for- 
mation in nearby galaxies and star-forming regions within 
galaxies rely on the Paa luminosity to interpret and to cal- 
ibrate othe^JjjFT^^djCjtor^^ 

(e.g.. ICalzetti et al.ll2005L l2007t lAlonso-Herrero et al.l l2006t 
iKennicutt et alj2007l;lDfaz-Santos et al.l2008l) . For extremely 
dusty, star-forming IR luminous galaxies with A(V) ~10- 
30 mag (e.g., Murphy et al. 2001; Dannerbauer et al. 2005; 
Armus et al. 2007), the extinction at Paa is <2 mag, and op- 
tically thin. Thus, even in the most obscured, luminous star- 
forming regions, the Paa line traces the intrinsic ionization 
rate, and thus the SFR. 

While a few studies have compared SFR indicators in dis- 
tant galaxies to the Ha emission line (see references above), 
no attempt to use Paa as a SFR indicator has yet been at- 
tempted in any galaxy of significant redshift. At z>2 the Paa 
A1.875 /im line shifts to A b s > 5.2 /im, and is accessible 
to the Infrared Spectrograph (IRS) on-board Spitzer, which 
covers the wavelength range of 5.2-38 /im. However, star- 
forming galaxies at z>2 are too faint intrinsically for obser- 
vations with IRS at the expected wavelength for Paa (5.4- 
8 /im), except for rare, extremely luminous objects, in which 
AGN likely dominate the H I line emission (e.g., Brand et al. 
2006), complicating any SFR-tracer comparisons. We have 
initiated a program using SpitzerflRS to study the Paa emis- 
sion in 12 galaxies at z > 2 gravitationally lensed by fore- 
ground galaxies or clusters of galaxies. Here, we discuss the 
detection of Paa in a star-forming galaxy at z = 2.515, grav- 
itatio nally lensed by the rich cluster Abell 2218 dKneib et al] 
120041) . This is the highest-redshift object with a measurement 
of the Paa emission line yet obtained, and thus we are able to 
study SFR indicators in a more "typical" star-forming galaxy 
directly at high redshift. 



The outline for the rest of this paper is as follows. In § 2, 
we summarize the properties of SMM J163554. 2+661225. In 
§ 3, we discuss the Spitzer observations and data reduction. In 
§ 4, we discuss the analysis of the spectroscopic and imaging 
data to derive rest-frame luminosities. In § 5, we compare the 
rest-frame luminosities of SMM J 163554.2+66 1225 to local 
samples, and discuss the utility of the various quantities as 
SFR indicators. In § 6, we present our conclusions. In the 
Appendix, we use stellar population models to analyze the 
rest-frame UV to near-IR spectral energy distribution (SED) 
of SMM J163554. 2+661225, and make some conclusions on 
the nature of the stellar populations in this galaxy. To derive 
physical quantities we use a cosmological model with Ho = 
70 km s" 1 Mpc" 1 , fi M , = 0.3, fi A ,o = 0.7. 

2. SMMT163554.2+661225 

We targeted SMM J 163554.2+66 1225 (16 /! 35 m 54.2 s , 
+66°12'24.5", J2000), which is a galaxy at z = 2.515 grav- 
itationally lensed by the galaxy clu ster Abell 2218. SMM 
J163554.2+661225 was identified bv lkneib et al.l J2004L here- 
after K04) as a lensed-galaxy candidate selected as a sub-mm 
source at 450 and 850 /im, with a very red optical counterpart. 
This source consists of a triply lensed system, with a com- 
bined magnification from gravitational lensing /i ~ 45 (K04). 
The brightest component (component "B" in the notation of 
K04) has /i = 22 ±2, K s = 19.5 mag with I-K s = 3.67 mag 
in the Vega-based system, and submillimeter flux densities 
of S450 = 75 ± 15 mJy and Ssso = 17 ± 2 mJy. The intrinsic 
sub-mm flux density (corrected for gravitational lensing mag- 
nification) is $850 = 0.77 mJy, a factor of two fainter than the 
sensitivity of the deepest sub-mm field surveys, and nearly an 
order of magnitude fainter tha n "typical" sub-mm-selected 
galaxies (IChapman et al.l 120051) . K04 provided the spectro- 
scopic redshift z = 2.5165 ±0.0015 based on the Ha emission 
line in a near-IR spectrum and z = 2.5 14 ± 0.001 from the ISM 
absorption lines in an optical spectrum (rest-frame UV). Fur- 
thermore, SMM J163554.2+661225 appears to be part of a 
larger group at z ~ 2.5, with at least two other galaxies lensed 
by Abell 2218 with z = 2.514 based on their UV absorption 
lines. 

iRigbv et alj d2008l hereafter R08) observed SMM 
J163554. 2+661225 with Spitzer as part of the time allo- 
cated to the guaranteed time observers (GTO). The GTO 
observations included imaging of Abell 2218 at 3.6-8 /im 
from the Infrared Array Camera (IRAC, iFazio et all 120041) . 
and imaging at 24-160 /im us ing the Multiband Imaging 
Photometry for Spitzer (MIPS, iRieke et all 120041) . In the 
Appendix, we use the available photometry from K04 and 
R08 to constrain the nature of the stellar populations and to 
estimate the stellar mass of SMM J163554.2+661225. R08 
targeted component "B" of SMM J1 63554.2+661225 for 
mid-IR spectroscopy with Spitzer/IRS dHouck et al.l l2004ab 
because this source has the brightest IR flux densities, with 
S24 = 1.16 ±0.10 mJy. Their IRS spectroscopy used the 
LL1 module, covering 19-38 /i m with R ~ 100 resolu tion. 
R08 fit templa t e IR SEDs from iDale & Heloul (120021) and 
ICharv & Elbazl (120011) to the 24, 450, and 850 /im data (in- 
cluding a MIPS 70 /im flux-density limit) to derive a total IR 
luminosity integrated from 8-1000 /im, Lir = L(8- 1000/im). 
They found L IR = (7.6 ± 1.9) x 10 11 L s , corrected for the 
gravitational-lensing magnification, where the range reflects 
differences in the template SEDs. This IR luminosity corre- 
sponds to a SFR, t/) = 140 ± 30 M Q yr" 1 , using the relation 
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FIG. 1.— Spitzer IRS observations of SMM J163554.2+661225. The 
grayscale image shows the IRAC 5.8 fim image of Abell 2218. The 
four IRS SL2 slit positions and orientation for the observations of SMM 
J163554. 2+661225 are indicated by the overlapping rectangles, centered on 
SMM J163554.2+661225 (component "B"). The horizontal black bar indi- 
cates a distance of 10 arcseconds. Arrows indicate component "B" (the IRS 
spectroscopic target) and a component "A", a fainter counter image of this 
galaxy (using the notation of K04). 

in iKennicuttl d 19981) for a Salpeter-like initial mass function 
(IMF) with lower and upper mass cutoffs of 0.1 and 100M Q , 
respectively. 

SMM J 163554.2+66 1225 appears to be undergoing an 
intense episode of star-formation with no evidence for 
an active galactic nucleus (AGN). K04 note that the ra- 
tio [N II]/Ha=0.3 ± 0.1 is consistent with ionization from 
star-formation for stars of near solar metallicity (see 
iPettini & Pagell 120041) . R08 show that the mid-IR spectrum 
contains strong emission attributed to polycyclic aromatic 
hydrocarbon (PAH) molecules, associated with dust heating 
from star-formation, with features similar to local starburst 
galaxies. The lack of any apparent silicon absorption at 

9.7 /j,m in the spec trum argues again st the presence of an ob- 
scured AGN (e.g., ISpoon et al.ll2007b . R08 analyzed archival 
data from the Chandra X-ray Observatory, which provided 
limits on the X-ray emission for this galaxy. This flux up- 
per limit corresponds to a luminosity limit of L(0.5-8keV) < 
5.0 x 10 42 erg s" 1 , corrected for the gravitational-lensing mag- 
nification. Compared to other sub-mm galaxies, the X-ray 
upper limit ru les out the presence of an AGN unless it is very 
obscured (e.g. lAlexander et ai1l2005l) . 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. IRS SL2 Spectra 

For our Spitzer/IRS observations, we targeted component 
"B" of SMM J163554.2+661225, using the IRS SL2 module, 
covering 5.2-8.7 /im. This is the same object targeted by R08, 
who observed this galaxy with the IRS LL module. Thus, the 
observations here extend the spectral coverage of this target to 
the shorter wavelengths, including the expected wavelength of 
Paa. 

The observed flux density of SMM J163554.2+661225 at 

5.8 pm is S 5 . 8 = 92 ± 9 pJy (R08). For the redshift of our tar- 
get, z = 2.515, we expect the Paa line to fall at 6.6 pm. We 
obtained our IRS observations in Mapping mode, placing the 
object at four different positions along the SL2 slit separated 
by 14". Figure Q] shows the position of the four SL2 slit po- 



sitions on the IRAC 5.8 fim image of the field for the epoch 
of observations. We observed SMM J163554.2+661225 us- 
ing 60 cycles of 60 s integrations, which we repeated at each 
of the four slit positions. The total integration time is 4 hours, 
for which the Spitzer Performance Estimation Tool predicts 
a signal-to-noise ratio of ~6-7 per resolution element in the 
continuum. In practice, we find that systematics reduce this 
slightly (see below). 

Observations were obtained on 2008 Feb 22-23 during a 
single AOR in order to prevent alignment uncertainties from 
different campaigns possibly at very different spacecraft ori- 
entations. We used high accuracy peak-up centroiding with 
the IRS blue filter to minimize any deviations in the absolute 
telescope pointing. We chose a peak-up star from the 2MASS 
catalog, against which we had matched the astrometric coor- 
dinates for SMM J163554.2+661225. 

We reduced the data starting with the S 17.2.0 Spitzer IRS 
pipeline data, which produced basic calibrated data (BCD) 
files. Our data contained some very strong cosmic rays 
which adversely affected the stray light correction in the IRS 
pipeline. However, we determined that the stray light effects 
are negligible in our data as the peak-up arrays do not con- 
tain bright sources during our exposures. We therefore used 
the BCD data products without the stray light correction (the 
BCD files f 2ap . fits and f 2unc . fits). 

The dominant background component at 5-8 pm is Zodia- 
cal light. Our observations included very long integrations of 
faint sources (these are currently the deepest SL2 observations 
taken of any object of which we are aware). Therefore we in- 
clude some additional steps to remove the sky background and 
detector effects. We generally follow the post-BCD reduc- 
tion steps used bv lTeplitz et al.l (12007b and lPope etaT] d2008l) 
in their analysis of deep IRS/LL data. We looked for indica- 
tions of latent charge in the array during our long exposures, 
as latent char ge accumulation is o bserved in long IRS/LL ex- 
posures (e.g.. iTeplitz et al.ll2007l) . However, we observe no 
increase in the median number of ADU s -1 over the course of 
a 60 min SL2 observation. Thus, there is no evidence for ap- 
preciable latent charge in the SL2 array over our integrations. 

We next identified and cleaned known bad and hot pixels 
using the SSC task IRSCLEAN 4 with the known warm-pixel 
mask for our IRS campaign (IRSX008900). We also identi- 
fied other "rogue" and "warm" pixels as pixels with abnor- 
mally high variance. We performed this latter task using an 
automated routine which computed the variance of each pixel 
over all the BCDs of an observation. We flagged rogue pixels 
as those pixels with high variance (> 4 x a 2 ). We then used 
IRSCLEAN to interpolate over these pixels. 

We constructed a sky spectrum for SMM 
J 163554.2+66 1225 at each of the four slit positions by 
combining the BCDs of the three other slit positions. To 
create the sky image, we took the median of the stack of 
each pixel after rejecting outliers using a sigma-clipping 
algorithm. We performed this process iteratively, masking out 
the location of SMM J163554.2+661225 during subsequent 
iterations (SMM J163554.2+661225 is the only source we 
identify in the 2D spectrum). We then subtracted the sky 
frame from each BCD and coadded the BCDs at each slit 
position. As a last step, we reran IRSCLEAN on the com- 
bined images for each slit position to clean any remaining hot 
pixels. These steps produced four 2D spectroscopic images 

4 http://ssc.spitzer.caltech.edu/postbcd/irsclean.html 
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FIG. 2. — Extracted Spitzer IRS SL2 spectrum of SMM J 163554.2+661225, covering 5.0-8.4 /im. The location of the Paa A 1.8751 emission line and possible 
CO absorption at a rest-frame of z = 2.515 are indicated. The dot-dashed line shows the continuum level predicted from the IRAC 5.8 and 7.9 pm flux densities, 
which is consistent with the measured continuum, excluding regions around the Paa emission line and excluding the possible CO absorption feature. The red- 
shaded spectrum shows the SL2 2nd order spectrum. The blue spectrum shows the SL2 1st order "bonus" spectrum. The error bars correspond to the uncertainties 
measured from the IRS observations, propagated through the data reduction, spectral extraction, and combination procedures. 



for SMM J 163554.2+66 1225, one at each of the four slit 
positions. 

We extracted ID spectra at each slit position using the 
Spitzer/IRS custom extraction (SPICE) software 5 with an op- 
timal extraction window. We extracted spectra from both the 
2nd spectral order (covering 5.2-7.7 /im) and the "bonus" 1st 
spectral order (covering 7.3-8.7 /im). SPICE provides the ex- 
tracted ID science spectrum and the propagated errors from 
the IRS observation. For each science spectrum, we also ex- 
tracted a sky spectrum using the same SPICE parameters, off- 
set from SMM J 163554. 2+66 1225. These sky spectra pro- 
vide consistency estimates for the error derived from the sci- 
ence spectra. We found these to be consistent, and we adopted 
the latter here in subsequent analysis. We then combined the 
four ID spectra as a weighted mean, using the weights de- 
rived from the error spectra. The estimate of the variance on 
the combined ID science spectrum is the inverse sum of the 
weights. Figure |2] shows the combined ID spectrum. 

SMM J163554.2+661225 is marginally resolved in the 
IRAC 5.8 and 8.0 fim images (channels 3 and 4). We es- 
timated the angular resolution using three stars in the field 
identified from HST/ACS imaging. For these three stars 
we measured a mean FWHM=1."8 in the channel 3 image. 
In comparison, for SMM J 163554.2+66 1225 we measure a 
FWHM=2."0 in this image. Most of the elongation of SMM 
J 163554.2+66 1225 appears to be aligned with the orienta- 
tion of the SL2 slit for our observations, as illustrated in fig- 

5 http://ssc.spitzer.caltech.edu/postbcd/spice.html 



ureQ] Although SMM J163554.2+661225 appears resolved, 
for the spectral extraction and flux calibration we neverthe- 
less assumed SMM J163554. 2+661225 to be a point source. 
This is reasonable because IRS SL2 has large pixels (1.8" 
pixels) compared to IRAC (1.2" pixels), so the spatial reso- 
lution is substantially worse compared to IRAC, and SMM 
J 163554.2+66 1225 is at best marginally resolved. In princi- 
ple, if it is resolved spatially, this could affect the flux cal- 
ibration of the extracted spectrum. However, we measure a 
continuum of SMM J 163554.2+66 1225 from our SL2 spec- 
trum that is within 1% of that expected from the measured 
IRAC 5.8 and 7.9 fim flux densities (see figure |2). Because 
pointing-induced slit throughput variations may produce un- 
certainties in the absolute flux measurements of 10-15%, we 
consider this agreement fortuitous. Nevertheless, we con- 
clude that our extracted spectrum corresponds to the spa- 
tially integrated emission of SMM J163554.2+661225, and 
we make no corrections because of the extended nature of this 
source. Furthermore, because the expected flux densities and 
measured continuum are essentially equal, we also make no 
additional corrections to the measured flux for light falling 
outside the SL2 slit. 

3.2. MIPS 70 fim Imaging 

SMM J 163554.2+66 1225 has SpitzerMlPS 70 and 160 fim 
data from the GTO observations. Using these, R08 placed an 
upper limit on the 70 fim emission of Sjq < 7 mJy. Since 
then, deeper MIPS 70 fim data have become available as part 
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of Spitzer program ID 30823. These new observations reach 
an exposure time of 5245 s at the deepest point, substantially 
longer than the GTO observations. We obtained the raw MIPS 
70 /im images from the Spitzer a rchive and reduced t hem us- 
ing the GTO Data Analy sis Tool dGordon et alj|2005l) follow- 
ing the steps described in lDole et al.l (120041) . We have applied 
further steps to mask sources detected in the image during the 
background subtraction steps, which greatly improves the fi- 
nal image quality. The final image achieves a limiting flux 
density of 070=0. 5 mJy for point sources. 

We detected point sources in the image using a weighted 
detection map and performed point-response function (PRE) 
weighted photometry. Because there are not enough bright 
point sources in the image to construct an empirical PRE 
from the image, we used the empirical PRF constructed 
from sources in wid er area (0.5 deg 2 ) extragalactic fields 
dPapovich et alj[2007b . Photometry is performed simultane- 
ously on c rowded s ource s using a version of the DAOPHOT 
software (IStetson| |1987i) . following the procedure from 
iPapovich et al.l (120071) . 

SMM J163554.2+661225 is marginally confused in the 
70 /im image. Figure [3] shows the contours of the 70 /im 
image overlaid on the 24 /im image. Given the substan- 
tially poorer resolution at 70 /im, the detected source cor- 
responding to SMM J 163554.2+66 1225 is likely blended 
with adjacent sources. We measured a flux density of Sjq = 
2.56 ± 0.90 for the source most closely associated with SMM 
J 163554.2+66 1225 (uncorrected for gravitational magnifica- 
tion). At the 70 /im angular resolution of Spitzer, SMM 
J 163554.2+66 1225 is blended with a counter image of itself 
(component "A" in the notation of K04), and possibly another 
source. Given that the ratio of the 24 /im flux densities of 
these components is S24(B)/S24(A)=1.6, we expect a similar 
ratio at 70 /im because presumably component A and B have 
the same 24 /im-to-70 fjm flux ratios (although K04 note that 
component B has a redder I-K color, and thus could have a 
higher S70/S24). To improve our flux-density measurements, 
we simultaneously measured point-source photometry in the 
70 /im image using a priori positions for sources detected in 
the 24 /im image. However, this method was very sensitive 
to uncertainties in the absolute astrometry between the 24 and 
70 /im images, and yielded a flux density for component B 
ranging from 670 = 0.0 mJy to 1.5 mJy (±1.0 mJy). There- 
fore in the analysis below, we use the direct photometry, and 
ascribe all the 70 /im emission to SMM J 163554.2+66 1225 
(component B). We note that the true 70 /im flux density may 
be lower by as much as 40% assuming the S70/S24 ratio above. 
As we will discuss below, the measured 70 /im flux den- 
sity implies a rest-frame L(24/im)-to-L(Paa) ratio lower than 
that for local IR-luminous galaxies with comparable bolomet- 
ric luminosities. Therefore, our assumption makes this result 
conservative. 

4. ANALYSIS 
4.1. Paschen-a Emission 

Figure 12 shows the combined IRS/SL2 spectrum. The dot- 
dashed line in the spectrum shows the predicted continuum 
from the measured IRAC photometry at 5.8 and 8.0 /im, S5.8 = 
92 /iJy and S7.9 = 110 /iJy, which agrees with the measured 
continuum to within 1%. The median signal-to-noise ratio of 
the continuum is S/N«5 from 5.5-7 /im. 

We identify the emission line at 6.591 /im as Paa at z = 
2.515, consistent with the expected redshift derived by K04. 
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FIG. 3.— SpitzerfMlPS 24 fim image of Abell 2218. SMM 
J163554. 2+661225 (Comp B, in the notation of K04) and its counter image 
(Comp A) are labeled. Contours show the Spitzer/MIPS 70 /im flux density, 
at levels of 1, 2, 3, 6 X the rms noise. Diamonds denote sources detected 
in the 70 fim data. The 70 (im source most closely associated with SMM 
J163554. 2+661225 is blended with its counter image and partially with the 
other nearby source. 

This is the highest redshift d etection of Paa in any g alaxy 
of which we are aware (cf., ISiana et al.l I2008L [2009). To 
measure the line parameters, we simultaneously fit the con- 
tinuum and the line, where we model the latter as a Gaus- 
sian. The Gaussian fit has a measured full-width at half max- 
imum, FWHM=0.067 ± 0.014 /im. This is comparable to the 
IRS/SL2 resolution, therefore we conclude the line is unre- 
solved. To derive the significance of this detection we gener- 
ated a series (10 3 ) of simulated spectra using the data and a 
random value taken from a Gaussian distribution with a equal 
to the derived uncertainty on each data point. We then remea- 
sure the line flux in each simulated spectrum, and take the 
inner 68% of the simulated distribution for each parameter as 
the uncertainty. The measured quantities are given in Table [TJ 
We derive a redshift for the Paa line of z = 2.515 ±0.003, 
which agrees within the measured uncertainties of K04 and 
R08 using spectral features in the optical, near-IR, and mid- 
IR wavelength ranges. 

We measure a Paa line flux of (8.6 ± 1.2) x 10~ 16 erg s -1 
cm" 2 , uncorrected for the gravitational lensing magnification 
or dust attenuation. We estimate the amount of dust attenu- 
ation affecting the nebular gas by comparing the ratio of the 
Paa line flux to the Ha line flux provided by K04 using K- 
band spectroscopy. Note that the K04 measurement of Ha has 
not been corrected for extended emission beyond their spec- 
troscopic slit of width 0"76. We estimate that the slit width 
used in the TT-band spectroscopy of K04 would miss up to 
35% of the light from SMM J163554.2+661225. If we ap- 
plied this (maximal) correction to the Ha flux, then it would 
decrease the extinction A(Paa) by 0.03 mag, and decrease the 
extinction-corrected Paa luminosity by «2%. Because these 
corrections are small, we do not apply any correction for the 
light fa lling outside the slit of the K-band spectroscopy. Us- 
ing the lCalzetti et al.l d2000l) attenuation law, we estimate the 
extinction to be A(V) = 3.6 ± 0.4 mag, which corresponds to 



6 



PAPOVICH ET AL. 



TABLE 1 
Derived Quantities 



A c EW F(Paa) L(Pact) A(V) A(Pact) L(Pact) col SFR 

(fim) z (A) (10~ 16 ergs" 1 cm 4 ) (lO^ergs" 1 ) (mag) (mag) (lO^ergs" 1 ) (M yr 1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 

6.591 ±0.006 2.515 ±0.003 363 ± 56 8.5 ±1.4 2.05 ±0.33 3.6 ±0.4 0.27 ± 0.03 2.57 ±0.43 171 ±28 



NOTE. — (1) Centroid wavelength of Pact line, (2) measured redshift, (3) rest-frame equivalent width of the Pact line, (4) measured line 
flux, with no correction for the gravitational lensing magnification or dust extinction, (5) lin e luminosity, correct ed for the gravitational lensing 
magnification fi = 22 (K04), (6) nebular extinction from Ha and Pact measurements using the Calzetti et al. 1 2000) law, (7) corresponding extinc- 
tion at 1.8751 fim, (8) line luminosity, corrected for the gravitational lensing magnification and dust extinction, (9) total SFR, corrected for the 
gravitational lensing magnification and dust extinction. 



A(Paa) = 0.27 ± 0.03 mag. Other attenuation laws that we 
tested (IRieke & Lebofskvll 19851: ICardelli, Clayton. & Mathisl 
[T989t iDraind 119891: iDopita et al.1 120051) correspond gener- 
ally to higher extinction est imates by as much a s 5A(V) w 
0.5 mag compared t o that of ICalzetti et al.l (120001) . However, 
ICalzetti et al.l J20 05) showed that compared with other atten- 
uation laws their attenuation law reproduces better the UV 
colors of starforming H II regions with similar extinction and 
properties as that derived here for SMM J163554.2+661225. 
Therefore, we adopt the Calzetti et al. law for our analysis 
here. As we will discuss below, the ratio of the mid-IR lumi- 
nosity to the extinction-corrected L(Paa) cor is lower relative 
to local IR galaxies of comparable luminosity. Using other 
attenuation laws would increase the Paa line luminosity by 
«30%. Therefore using the lower extinction value provided 
by the Calzetti law also provides a conservative choice. Nev- 
ertheless, the dust-extinction correction remains a systematic 
uncertainty. 6 

We derive a Paa line luminosity of L(Paa) cor = (2.57 ± 
0.43) x 10 42 erg s" 1 after applying the extinction correction, 
2()0.4A(Paa) = 1.28, the gravitational lensing magnification, fi = 
22, and the luminosity distance, di = 2.054 x 10 4 Mpc, at z = 
2.515 for the default cosmology. This luminosity corresponds 
to an ionizing continuum flux of Q(H°) = (1 .6 ± 0.3) x 10 55 7 
s _1 (lQsterbrocklll989t lKennicuttl[l998l) . For an IMF with a 
Salpeter-like slope from 0.1 to 100 Mq, the implied SFR is 
171 ±28 Mq yr" 1 . The statistical uncertainty on the SFR is 
17%, which is the highest accuracy on a SFR derived for a 
high-redshift galaxy to date. Note that while this uncertainty 
neglects the uncertainty on the gravitational lensing magnifi- 
cation (w 10%, see § 2), the lensing is expected to be achro- 
matic, and it will systematically scale all derived luminosities 
and SFRs. 

4.2. CO Absorption 

The IRS/SL2 lst-order "bonus" spectrum shows a possible 
absorption feature at 8-8.5 /im(see figure©. Atz = 2. 515 this 
corresponds to «2. 3-2.4 fim, and would imply strong molec- 
ular CO absorption. If this feature is real, it would be the 
first time it has been observed in a galaxy with any signifi- 
cant redshift. CO absorption occurs in the atmospheres of red 
supergiant stars (K- and M-type), primarily from post-main- 
sequence O-type stars (e.g., Doyon, Joseph, & Wright 1994; 
Ridgway, Wynn- Williams, & Becklin 1994; Goldader et al. 

6 We have learned that a recent reanalysis of the K04 near-IR spectroscopy 
(Richard et al., in preparation) provides a Hct flux that is Ri4x fainter than 
that of K04. If this flux is correct, it would increase the extinction to A(Hct) = 
4.9 mag and A(Pact) = 0.45 mag, and would increase the extinction-corrected 
Pact luminosity (and SFR) by « 18%. However, none of the main conclusions 
here would be affected significantly. 



1995). Given the quality of our SL2 spectrum, a detailed anal- 
ysis of the CO absorption is cautionary. Nevertheless, we de- 
rive a spectr oscopic CO index of CO sp = . 29 ± . 05 using the 
definition of iDoyon et alJ (11994 see also Smith et al. 1996). 
Combined with the slope of the rest-frame near-IR contin- 
uum, this places SMM J163554.2+661225 in the "[Dust- 
1 Red dened Starburst Population" locus of iRidgwav et alJ 
(119941) . implying that the starburst in this galaxy dominates 
the rest-frame near-IR emission. 

If the strength of the CO absorption feature in SMM 
J163554. 2+661225 is real, then the age of the starburst cor- 
responds to the time it takes O-type stars to enter the red su- 
pergiant phase. The fact that the H I recombination lines re- 
main strong requires ongoing early-type (primarily O-type) 
star formation. The combination of these facts impli es an age 
for the starburst of ^10-50 Myr (IDoyon et alJll994l) . 

4.3. The Mid-Infrared Luminosity 

Using the MIPS 24 and 70 fim, and IRS/LL observations, 
we derive monochromatic luminosities at rest-frame 8 and 
24 fim, L(8/im) = ^L„(8/im) and L(24^m) = ^L„(24/im). 
Below we compare our results for SMM J163554.2+661225 
against datasets of local galaxies where the above quantities 
are derived using photometry measured from IRAC 8 fim and 
MIPS 24 fim. Therefore, we compute L(8/im) and L(24fim) 
which match the rest-frame IRAC and MIPS bandpasses as 
closely as possible. 

At z = 2.515 the MIPS 24 fim photometry corresponds to 
rest-frame 6.8 fim. To convert this to rest-frame IRAC 8 fim 
we use the IRS/LL spectrum of SMM J 163554.2+66 1225 
from 20-38 fim from R08, corresponding to rest-frame 5.7- 
10.8 fim. Figure |] shows the IRS/LL spectrum of SMM 
J 163554.2+66 1225 from R08 overlaid with the MIPS 24 fim 
bandpass in the observed frame and the IRAC 8 fim bandpass 
in the rest frame. We integrate the IRAC 8 fim transmission 
function with the IRS/LL spectrum (shifted to the rest-frame) 
to derive a rest-frame 8 fim flux density, /„ o(8/^m)=1.5 ± 
0.2 mJy, where the subscript "0" denotes the rest-frame quan- 
tity. This corresponds to /„ o(8/^m)=70± 10 fiJy corrected for 
the gravitational magnification. In contrast, R08 measured a 
24 fim flux density of 524 = 1-16 mJy. This difference is pri- 
marily a consequence of the fact that at z = 2.515 the MIPS 
24 fim bandpass is mostly insensitive to the strong 7.7 fim 
and 8.6 fim PAH features, which dominate the mid-IR spec- 
trum of SMM J 163554.2+66 1225. These PAH features are 
included in the rest-frame IRAC 8 fim bandpass, which re- 
sults in the much higher flux density when averaged over this 
bandpass. 

We measured the total IR luminosity for SMM 
J 163554.2+66 1225 by fitting a suite of template IR SEDs 
dCharv & Elbazll200U iDale & Heloull2002t IRieke et al.ll2009l) 
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to the 70, 450, and 850 fim photometry, as illustrated in 
figure [5] (see also § 5.1). The best-fit templates corre- 
spond to a total IR luminosity, Lir = (5 - 10) x 10 11 Lq, 
corrected for the gravitational lensing magnification. We 
estimated the rest-frame 24 //m flux density for SMM 
J163554.2+661225 by integrating these best-fit IR SEDs 
with the MIPS 24 /im transmission function. This yielded, 
/„ i o(24/im) = 0.24 ± 0.09 mJy, corrected for the gravitational 
lensing magnification. There is an additional small systematic 
uncertainty resulting from differences in the template IR 
SEDs, corresponding to cr sys = 0.03 mJy. The error bar here 
is dominated by the uncertainty on the 70 fjm photometry, 
which is the closest band to the rest-frame 24 /im datum. 
Furthermore, we remind the reader that the 70 / im flux 
density may be lower by as much as 40% (see § 13.2b . which 
would decrease the 24 /im luminosity. 

We convert the rest-frame flux densities to luminosities us- 
ing, L(A) = ^/i/,o(A) x {A-Tid}) x fjT 1 , where /t is the gravi- 
tational lensing magnification, and cIl is the luminosity dis- 
tance. Applying this formula to the 8 and 24 /im rest- 
frame quantities derived above yields L(8/mi) = (3.8 ±0.6) x 
10 44 erg s" 1 and L(24/im) = (4.3 ± 1.6) x 10 44 erg s" 1 . 

5. DISCUSSION 

The observations for SMM J163554.2+661225 provide in- 
dependent estimates for the total SFR. In particular, we com- 
pare the Paa luminosity which stems from the ionized gas 
in H II regions and traces the number of ionizing photons, to 
estimates from the mid-IR and total IR luminosities, which 
measure prima rily the dust-reprocessed emission fr om mas- 
sive stars (e.g.. lKennicun1ll998trKewley et al.ll2002l) . As dis- 
cussed in § 4. 1, the Pact line luminosity corresponds to a SFR 
^=171±28M Q yr 1 . 

5.1. The Pact Luminosity Compared to the Total IR 
Luminosity 

Local star-forming galaxies and star-forming regions 
show a tight correlation between L(Paa) and Lir (e.g., 
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FIG. 4.— IRS/LL spectrum from 20-35 ^m for SMM J163554.2+661225 
from R08. The IRS spectrum is shown as the yellow-filled histogram (R08). 
The blue-hashed region denotes the MIPS 24 fim bandpass in the observed 
frame at z = 2.515. The red-hashed region denotes the IRAC 8 (im bandpass 
in the rest-frame. We integrated the IRS/LL spectrum with the 8 (im band- 
pass in the rest-frame to derive the rest-frame 8 fim flux density. This differs 
significantly from the observed MIPS 24 (im flux density because of the con- 
tribution from strong spectral features in this wavelength region attributed to 
PAHs, especially at 6.2, 7.7 and 8.6 ^m. 
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FIG. 5.— Infrared SED of SMM J163554.2+661225. In the top panel, the 
data points show the MIPS 24 and 70 /im, and SCUBA 450 and 850 /im mea- 
sured flux densities. The right axis shows the flux densities after correcting 
for the gravitational lensing magnification. The top axis shows the rest-frame 
wavelength for z = 2.515. The curves show IR SE D template fi t s to th e 70, 
450, and 85 flux densi t ies, us ing templates from Rieke et al. 1 20 0^, solid 
black line), IDale & Heloul <250l dashed line), and Charv & Elbai OOP"! 
dot-dashed line). Formally the Chary & Elbaz and Rieke et al. templates 
provide better fits to the data. However, uncertainty in the total IR luminosity 
on the order of a factor of two remains owing to differences in the templates. 
The range of the implied total IR luminosity, Lir = L(8- 1000/im), ranges 
from 5- 10 X 10" Lq, depending on the model adopted. To improve these 
constraints requires flux density measurements at observed wavelengths of 
~ 150-250 /an, which will constrain the peak of the thermal dust emission. 
The bottom panel shows the ratio of the best-fit models and data points to the 
best-fit model of Dale & Helou. Note that all of the best-fit models to the 
far-IR data underpredict the observed 24 (im flux density by factors of >2. 



ICalzetti et al.ll200H lAlonso-Herrero et al.ll2006h . For SMM 
J 163554.2+66 1225 R08 derived a total IR luminosity in the 
range Lir = 5.7-9.5 x 10 11 Lq, taking into account the sys- 
tematic uncertaint y owing to the choice of IR SEDs (see also 
the discussion in iPapovich et alj|2007h . Here, we reanalyze 
the total IR luminosity of SMM J163554.2+661225 by fitting 
different sets of IR SEDs to the flux density at 70, 450, and 
850 /im flux densities, all of which sample the thermal IR 
emission (see figure|5]and § 4.3). We exclude the 24 /im flux 
density from this analysis because it probes ~6-7 /im in the 
rest-frame mid-IR, and its relationship to the thermal far- 
IR emission is not straightforward. Figure [5] shows the best 
fit IR SED te mplates and measured flux densities. Using the 
IDale & Heloul (120021) templates, we derive L IR = (5.0 ±0.6) x 
10 11 Lq, with a goodness-of-fit, x~ jv = 3.2 an d v = 2. We find 
forma lly better fits using both the IR SEDs of Ch arv & Elbazl 
(jO p), giving Lir = ( 5.9 ± 0.6) x 10 11 L with X 2 /v = 0.8, 
and lRieke et all d2009t) . giving Lir = ( 1 ± 1 . 1 ) x 1 1 1 L with 
X 2 jv = 0.2. The IR luminosities from the fits are consistent 
with the results from R08, where the difference here is that we 
have excluded the 24 /im flux density and include the deeper 
70 /im flux-density measurement. Because the different IR 
SED templates all are consistent with the data, they imply 
there is a factor of 2 uncertainty on the total IR luminosity 
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FIG. 6. — Paa luminosity versus the monochromatic luminosity at rest-frame 8 fim. The left panel shows L(Pace) versus L(8/im) and the right panel shows 
L(Paa) versus the ratio L(8/im)/L(Paa). In both panels the symbols and lines are the same. The large, red-filled circle shows the measured value for SMM 
J163554. 2+661225. The light-red-shaded area indicates the error bar on the ratio for SMM J163554. 2+661225. Diamonds show local luminous IR galaxies 
(Lir = 10" -10 12 Lq) from the sample of lDi'az-Santos et al.1 <2008l DS08) Triangles show individual star-form ing regions in M51 ICalzetti et al.1 <2005l C05). 
Stars show starburst and low-metallicity galaxies from Calzetti et al. (2007, C07) and Enselbracht et al. 12005, E05). The dot-dashed line shows the best-fit 
relation to individual star-forming H II r egions in M51 from Calzetti et al. (2005|). The short-dashed line shows the best-fit r elation derived for star-fo rming 
galaxies and individual H II regions from Calzetti et al. 1 2007), which is similar the relation for luminous IR galaxies derived by Dfaz-Santos et al. ( 2008, long- 
dashed line). The 8 (im and Pace luminosities for SMM J163554. 2+661225 are consistent with the extrapolated relationship for local star-forming regions and 
star-forming galaxies. 



owing to differences in the choice of template. Moreover, be- 
cause the IR SEDs we have tested are not continuous in Lir, 
this implies there is at least a factor of 2 uncertainty on total 
IR luminosities of z ~ 2.5 galaxies, even when flux densities 
are available at 70 /im and sub-mm wavelengths. 

The total IR luminosities correspond to a range of SFR, 
tp = 90- 180 Mq yr. The upper range of the SFR corresponds 
to the fit using the Rieke et al. template, and these are con- 
sistent with the SFR derived from the Paa luminosity. The 
templates from Dale & Helou and Chary & Elbaz yield lower 
Lir values, and imply SFRs lower by ^40-50% compared to 
that from L(Paa). To improve the accuracy on the IR luminos- 
ity will require flux density measurements at observed wave- 
lengths of « 150-250 /im, in order to constrain the peak of the 
thermal dust emission (see figure[5]l. Nevertheless, the current 
analysis provides evidence that the total IR luminosity and 
Paa luminosities are consistent for SMM J163554.2+661225 
within a factor of 2. 

Parenthetically, we note that none of the empirical IR tem- 
plates are able simultaneously to fit both the thermal IR emis- 
sion and the strength of the PAH emission features in the 
mid-IR. This is apparent in the lower panel of figure [5] where 
the best-fit IR templates imply lower observed 24 /im flux 
densities less than the measured value. This issue has been 
identified in the anal ysis of other high- redshift gravitationally 
lensed galaxies (e.g. JSiana et al .1120081) . 

5.2. The Paa Luminosity Compared to the 8 fim Luminosity 

The origin of the mid-IR emission (rest-frame 5-10 /im) 
in local galaxies is attributed both to very small-grain dust 
continuum emission and molecular PAH emission. Both 
the PAHs and very small grains are heated both by ion- 
izing and non-ionizing sources (in particular, the am- 



bient galactic radiation field from evolved stars, Li & 
Draine 2002). Local galaxies and star-forming regions 
show a nonlinear relation between L(8/zm) and L(Paa ) 
(see, e.g.. iForster Schreiber et aT]|2004t ICalzetti et al.ll2007l) . 
where the slope of the correlation depends on the dust- 
heating source(s ), the dust spallation/formation rates, and the 
metallicity (e gjHouck et al.ll2004btlEngelbracht et aT] |200l 
ICalzetti et atll2007HDraine et al.ll2007l) . 

The MIPS 24 /im flux density probes rest-frame mid-IR, 
~6-7 /zm at z=2.515, as illustrated in figure|4] Many studies 
of the IR emission in distant galaxies use the measured 24 /im 
flux density to estimate the total IR emission, and we test these 
relations for SMM J163554.2+ 661225. Using the 24 / mi flux 
density with the prescription of iPapovich et al.l (|2006) yields 
a estimated total IR luminosity, L IR = (1.2 ±0.1) x 10 12 L©, 
where the error is statistical only and does not include sys- 
tematic uncertainties (see discussion in Papovich et al. 2006). 
This corresponds to a SFR, ip ~ 220 Mq yr" 1 . While this is 
consistent with the SFRs derived from the Paa emission and 
Lir measured from the far-IR data, this is somewhat fortu- 
itous because the IR template used to extrapolate the observed 
24 /zm flux density implies higher flux densities at observed- 
frame 70 and 850 /im compared to the observations. This is 
similar to the statement in §5.1 that none of the template IR 
SEDs are capable to fit simultaneously the far-IR flux densi- 
ties a nd mid-IR emission fe atures. Using the scaling relation 
from IPapovich et al.l (120071) . which includes bolometric cor- 
rections using the average 70 and 160 /tm flux densities of 
1.5 < z < 2.5 galaxies, yields a nearly equal estimate for Lir 
as using the 24 /im data only, with a similar offset compared 
to the implied SFR from Paa. 

R08 derive a scaling relation between the rest-frame 
8 /im luminosity, L(8/im), and Lir for their sample of 
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FIG. 7. — Pact luminosity versus the monochromatic luminosity at rest-frame 24 fim. The left panel shows L(Paa) versus L(24^tm) and the right panel 
shows L(Paa) versus the ratio L(24/xm)/L(Pact). In both panels the symbols and lines are the same. The measured value for SMM J163554. 2+661225 is 
shown as the large, red-filled circle. The light-red-shaded area indicates the error bar on the ratio for SMM J163554. 2+661225. Note that L(24^m) for SMM 
J163 554. 2+661225 may be low er by Ri40% owing to crowded photometry (see § 3.2). Squares show local ultraluminous IR galaxies (ULIRGs, Lir > 10 12 Lq) 
from ITTannerbauer et al. 12005, D05), with L(24/im) estimated using the method discussed in the text. Diamonds show local luminous IR galaxies (LIRGs, 
Lir = 10" - 10 12 Lq) from the sample offA] onso-Herrero et al. 12006. AH06); the diamond indicated by the circle corresponds to the galaxy IC 860, which is 
proble matic according to A06. The range of the ordinate in the right panel excludes IC 86 0. Triangl es show local star-formin g H II regions in M51 Calz etti et all 
120051 C05). Stars show starburst and low-metallicity galaxi es fromlCalzetti et a l. 12007, C07) and Engelbracht et al. l2005| E05). The short-dashed line shows 
the best fit relationsh ip to local luminous IR galaxies from Alonso-Herrero et al. 120061) . The long-dashed line shows the best-fit to local g alaxies and star- 
forming regions from Calzetti et al. 12007). The short-dashed line shows the best-fit relation for individual H II regions Calzetti et al. ('200$. The rest-frame 
24 fj,m luminosity for SMM J163554. 2+661225 is significantly lower (by ks 0.3-0.5 dex) with respect to local galaxies of similar Paa luminosity. In contrast, 
the L(24/xm)/L(Paa) for target is consistent with that found in individual star-forming regions. 



gravitationally lensed z ~ 2 galaxies, which includes SMM 
J163554.2+661225. Using the L(8/zm) derived in § 4.3, the 
R08 scaling relation yields Lir = (7.9 ± 1.5) x 10 11 L . This 
is consistent within the range of the total IR luminosity de- 
rived above, and consistent with the implied SFR from Paa 
w ithin the errors. 

iPope et all ( 120081) derive scaling relations between Lir 
and the luminosity of the 6.2 and 7.7 /xm PAH emission 
features. R08 measured their PAH luminosities for SMM 
J 163554. 2+66 1225 simultaneously using PAHFIT (Smith et 
al. 2007). However, as discussed by Pope et al., line fluxes 
from measured by PAHFIT are generally higher than those us- 
ing their method (see also the discussion in Sajina et al. 2007; 
Siana et al. 2008). Thus, we remeasured the 6.2 and 7.7 /m 
PAH luminosities individually, fitting each emission line with 
a Drude profile while fitting the slope and intercept of the 
continuum. Our fits for each line yielded L(6.2/im,PAH) = 

and L(7.7 y um,PAH) = (8.89 ± 
Using equations 4 and 5 from Pope et 
al. we infer Lir « 2.3 x 10 12 L Q and 3.1 x 10 12 L o for the 6.2 
and 7.7 /im PAH features, respectively. While these estimates 
of Lir are higher by a factor of order three compared to that 
measured from the far-IR data, they are within the scatter ob- 
served in the the relation between the PAH luminosities and 
total IR luminosity in the Pope et al. sample. However, the 
intrinsic IR luminosity of SMM J 163554. 2+66 1225 is a fac- 
tor of two lower than the high-redshift sub-mm from Pope et 
al., and it is possible the extrapolated rela tion does not hold . 
This may be supported by the results of ISiana et alJ (120081 
120091) . who observe a similar offset between the PAH lumi- 



nosity and total IR luminosity in their study of intrinsically 
less-luminous, lensed UV-bright objects. 

Figure [6] shows the Pace luminosity against the rest-frame 
8 yum luminosity for SMM J163554.2+661225, compared 
against luminosi ties for samples of local galaxies and star- 
forming regions (ICalzetti et al.ll2005Ll2007HEngelbracht etaD 



(2.25 ± 0.08) x 10 43 erg s" 1 
0.04) x 10 43 erg s" 1 . 



12005b iDfaz-Santos et alJl2008l) . ICalzetti et al] (120071) derive 
a scaling relation between the Paa and 8 /^m luminosi- 
ties, L(8/xm) oc L(Paa) Q , with a = 0.94. Note that Calzetti 
et al. derived this correlation in terms of luminosity sur- 
face densities (luminosit y per unit physical area). However, 
IDfaz-Santos et al] (120081) obtain a similar slope for the corre- 
lation between the 8 /im and Paa luminosities. While the lu- 
minosities for SMM J 163554.2+66 1225 appear broadly con- 
sistent with the extrapolated relations, a clearer picture is evi- 
dent by comparing the ratio of L(8yum)/L(Paa), shown in the 
right panel of Fig [6] There is much scatter in the local sam- 
ples, but L(8/rni)/L(Paa) for SMM J 163554.2+66 1225 agrees 
broadly with those ext rapolated relationships that include the 
star-f orming galaxies (ICalzetti et al. 200"7t IDfaz-Santos et al] 
12008b . Using only the extrapolated relation from individual 
H II regi ons would underpre dict the amount of 8 /im emis- 
sion (cf., ICalz etti et al. 2005). This implies that the fraction 
of photons from star-formation reradiated as 8 /im luminosity 
is weakly dependent on ionizing luminosity. 

5.3. The Pact Luminosity Compared to the 24 [im Luminosity 

The emission at rest-frame 24 /im results from thermal dust 
grains heated by ionizing and non-ionizing sources. Em- 
pirically, local galaxies and star-forming regions follow a 
correlation with L(24/im) oc L(Paa)", with a in the range 
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1.03-1. 23 (lCalzettietalJl200l 120071: lAlonso-Herrero et alJ 
120061) . lAlonso-Herrero et al.l (120061) argue that the correla- 
tion with a > 1 arises as dust absorbs ionizing and UV- 
continuum photons with increased efficiency in more heav- 
ily obscured, more luminous systems, so that an increas- 
ing fraction of the bolometric luminosity associated with 
star for mation emerges in the IR with warmer dust temper- 
atures JLonsdale Persson & Heloul 119871 : IWang & Heckmanl 
[l996t IDraine & Lil 120071) . ~ similar conclusion is reached 
by lCalzetti et all (120071) . who argued that the a > 1 correla- 
tion exists because objects with increasingly higher starlight 
intensity and L(Paa) have higher dust temperatures, where 
the peak of t h e emis sion moves to sho rt er IR wavelengths. 
ICalzetti et ail ( 12007b and IDraine & Lil d2007l) discuss the 
physical basis for this correlation. 

Figure [7] shows the Pact luminosity plotted against the 
rest-frame 24 fim luminosity derived abo ve for SMM 
J 163554.2+66 1225 and for local sam ples (ICalzetti et all 



2005 



2005 



120071: lEngefbracht et ail 120051: iDannerbauer et al.l 
lAlonso-Herrero et al.l 120061)" Interestingly, SMM 
J163554. 2+661225 has lower L(24/im) than local galaxies of 
comparable L(Pact), where the offsetis ps0.3-0.5 dex. This is 
more apparent when comparing the ratio of L(24/im)/L(Pact), 
shown in the right panel of figure [7] Furthermore, because 
the 70 fim data for SMM J 163554.2+66 1225 is blended (see 
§ 3.2), the rest-frame 24 /im luminosity may be lower than 
indicated in the figure, making the offset between SMM 
J 163554.2+66 1225 from the relation for the local samples 
more pronounced. 

SMM J163554.2+661225 has an estimated total IR lumi- 
nosity, Lir = (5- 10) x 10 11 Lq, comparable to local ultralumi- 
nous IR galaxies (ULIRGs, Lir > 10 12 Lq) . Figure |7Jincludes 
data fo r local ULIRGs from the sample of IDannerbauer et all 
(120051) . where we combine their Pact measurements with 
IRAS 25 fim measurements from the literature 7 . The Danner- 
bauer et al. measurements of Pact come from longslit near- 
IR spectroscopy, and we have made no attempt to correct for 
emission outside the slit. The Pact-emission in many of local 
ULIRGs likely results from very compact, nuclear regions, 
and therefore the spectroscopic slit should contain most of this 
emission. Nevertheless, we removed from the local ULIRG 
sample those objects with 2MASS isophotal diameters >15", 
and we also removed those objects with spectroscopic signa- 
tures of AGN. These steps excluded roughly one-third of the 
sample, including the most egregious outliers. Nevertheless 
we caution that significant uncertainty may remain due to pri- 
marily unknown Pact emission outside the spectroscopic slit 
or other aperture effects. 

Even with these caveats the local ULIRGs follow the ex- 
trapolation of the local L(24/im) and L(Pact) relation, but they 
show a large scatter. We suspect that the large scatter results 
for the reasons discussed above, and there may be additional 
components to the dust heating beyond ionization from early- 
type stars, including AGN and the ambient galactic emis- 
sion. However, only one galaxy in the local ULIRG sam- 
ple has L(24^m)/L(Pact) and L(Pact) comparable to SMM 
J 163554.2+66 1225, and this galax y (IRAS 04384-4848) ha s 
a highly uncertain dust correction (IDannerbauer et ai]l2005l) . 
This implies that no (or at best, few) low redshift ULIRGs 
have similar physical conditions producing comparable ratios 
of mid-IR-to-Pact luminosity. 

7 see Moshir, Kopman, & Conrow (1992); the NASA Extragalactic 
Database (NED), http://nedwww.ipac.caltech.edu/ 



100 



i 10 



c 

CD 
Q 



Rest-frame Wavelength [am] 

10 100 1000 

, ■ u 10 





A(V) = 2.2 ma. 
A(V) = 4.5 mag 
A(V) = 6.7 mag ^ 
A(V) = 9 mag 
A(V) = 18 mag ij, 

A(V) = 36 ma< 

A(V) = 72 ma- 




10 100 1000 

Observed Wavelength [fim] 



FIG. 8.— Infrar ed SEP of SMM J163554 . 2+661 225 compared to radiative 
transfer models of Siebenmorgen & Kriigel 1 2007). The data points show the 
measured flux densities at IRAC 3.6, 4.5, 5.8 and 8.0 /im, the MIPS 24 and 
70 fim, and the SCUBA 450 and 850 fim. The right axis shows the flux 
densities after correcting for the gravitational lensing magnification. The top 
axis shows the rest-frame wavelength forz = 2.515. The curves show model 
fits to the 24, 70, 450, and 850 fim flux densities. These SEDs are computed 
for a spherical PDR ionized uniformly by an interior starburst. The resulting 
amount of visual extinction is a variable in the model. Each curve shows the 
best-fit model as a function of visual extinction, as indicated in the figure 
inset. 



iKennicutt et al] (120091) combine Ha emission-line mea- 
surements (uncorrected for dust extinction) and IR contin- 
uum measurements of local star-forming galaxies, and derive 
SFR calibrations of the form, tp = 7.9 x 10" 42 x [L(Ha) obs + 
a\L(X)]. For the IRAC 8 and MIPS 24 fim rest-frame 
bands they obtain «g = 0.01 1 and «24 = 0.020. Using the mid- 
IR luminosities for SMM J16 3554.2+661225 der ived above, 
and L(Hct) obs = 3.7 x 10 8 L Q dKneib et all 12004 . we obtain 
ip ~ 45 and 80 M Q yr -1 , for the 8 fim and 24 fim luminosi- 
ties, respectively. These are lower by factors of 4 and 2 com- 
pared to that derived from the Pact luminosity, but the y are 
within the dispersion reported bv lKennicutt et alJ ( 120091) . The 
intrinsic SFR of SMM J 163554.2+66 1225 is also consider- 
ably larger than the objects used to calibrate these relations 
(see also Moustakas & Kennicutt 2006), and it is possible the 
calibrations do not apply under extrapolation. Larger samples 
of luminous, high-redshift galaxies are needed to test these 
relations. 



5.4. The Nature of SMM J '163554.2+661225 

SMM J163554.2+661225 appears to host heavily obscured 
star-formation at a rate, t/j w 170 M yr" 1 . The hydrogen ion- 
ization rate is go = (1.6 ±0.3) x 10 55 7 s" 1 , implying SMM 
J 163554.2+66 1225 m ay contain as many of ~ 10 6 O stars 
dSternberg et al.ll2003l) . Given the estimate for the molecular- 
gas mass (M[H 2 ] « 4.5 x 10 9 M Q , Kneib et al. 2005), this 
galaxy could sustain this SFR for t ~ 30 Myr. The starburst 
in SMM J 163554.2+66 1225 has had a duration of < 100 Myr 
based on the analysis of rest-frame UV-to-near-IR SED (see 
the Appendix) and supported by the strength of the possi- 
ble CO index (§ 4.2). This is comparable to the dynamical 
and gas-consumption timescales based on the observations of 
the molecular gas. Because the gas-consumption timescale 
is consistent with the dynamical time and the starburst age, 
SMM J 163554. 2+66 1225 is likely about midway through this 
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stage of enhanced star formation. 

The Pact-derived SFR for SMM J 163554.2+66 1225 is 
consistent with the SFR implied by the total IR luminos- 
ity and the rest-frame 8 [im luminosity. However, the 
rest-frame 24 p,m luminosity is significantly lower than 
that expected from the Pact luminosity based on the lo- 
cal relations. This implies that SMM J163554.2+661225 
lacks a warm (~100 K) thermal dust component typical 
of local /flAS-selected st ar-forming galaxies of compara- 
ble bolometric lumi nosity (lLonsdale Persson & Heloulll987l : 
ICalzetti et al.l 120001) . which drives the non-linear relation- 
ship betw een the mid-IR lum inosity and the Pact luminos- 
ity (e.g., ICalzetti et ail 120071) . Indeed, the thermal dust 
emission of SMM J163554.2+661225 (see figure |3) peaks 
near ~100 /xm, and is consistent with nearly pure emis- 
sion from dust with temperature, 7b = 52 K, and emissiv- 
ity, (3 = 1.5. This is very similar to the dust emission of 
photodissoci ation regions (PDRs) in the vicin i ty of H II re- 
gions (e.g.. lLonsdale Persson & Heloul 119871 : ICalzetti et alj 
I2000t IChurchwellll2002l) . In contrast, the IR SEDs for local 
galaxies with L IR comparable to SMM J163554.2+661225 
have significant contri butions of warm (> 70 K) dust to 
the IR emission (e.g.. ICharv & Elbazl I200U iDale & Heloul 
l2"002USiebenmorgen & Krugelll2007HRieke et al.ll2009l) . This 
is a qualitatively different than what is observed in SMM 
J163554. 2+661225. Indeed, the fact that no galaxies in the 
local ULIRG sample have L(24^tm)/L(Paa) values as low 
as SMM J163554.2+661225 implies that star formation in 
some high-redshift galaxies with ULIRG luminosities is fun- 
damentally different. The lower ratio, L(24/xm)/LiR ~ 0.1, 
for SMM J163554.2+661225 implies a lower starlight in- 
tensity than th at in local galaxies of comparable bolomet- 
ric luminosity dDraine & Lil 120071) . Therefore, while SMM 
J 163554.2+66 1225 appears to host a massive starburst that 
is similar to local star-forming regions, it is dramatically 
"scaled-up" in luminosity (and presumably SFR). 

This conclusion is supported by a comparison of the mid- 
IR and fa r-IR emission of SMM J16355 4.2+661225 to the 
models of Si ebenmorgen & Kr tigel (2007]), who calculate the 
IR emission for PDR-like regions with spherical symmetry of 
variable size surrounding H II regions ionized by a starburst 
with a variable fraction of the luminosity coming from OB 
associations. As illustrated in figure[8] the best-fitting models 
from Siebenmorgen & Kriigel correspond to an obscured stel- 
lar population where 60 - 90% of the luminosity originates 
from OB associations. Such models require emission from 
dust clouds with a range of extinction, A{V) 2-72 mag. 
Models with visual extinction A(V) < 18 mag that reproduce 
the data have an intrinsic luminosity ss 10 10 - 10 11 Lq, and 
this must be scaled up by a factors of wlO-100 to match 
the data, producing total IR luminosities in agreement with 
that derived above (§ 5.1). Scaling the size of the PDR 
according, this implies a radius for the ionization-front of 
~3-4 kpc, consistent with the spatial extent of the molec- 
ular CO emission of SMM J163554.2+661225 (~3 kpc x 
1.5 kpc. iKneib et al.ll2005l) . Models with higher intrinsic IR 
luminosity (;§> 10 11 Lq) do not reproduce the data. Mod- 
els with higher extinction require larger intrinsic luminosi- 
ties (as much as an order of magnitude larger for the model 
with A(V) = 72 mag), and these provide worse fits to the 
data (thus, they are less physical). Furthermore, models 
with A(V) > 30 mag would be optically thick to Pact pho- 
tons, and cause significant attenuation of the mid-IR emis- 



sion. Whi le such models ar e physically motivated in some 
cases (e.g.. lRiekeetaLll2009l) . if this were the case for SMM 
J 163554.2+66 1225 then the attenuation correction to the Pact 
luminosity would imply a much larger SFR, which should 
be substantiated by a larger total IR luminosity. The general 
agreement between the dust-corrected Pact luminosity and to- 
tal IR luminosity (§5.1) excludes models with very high ex- 
tinction. 

While the models assume a spherical shell-like dust config- 
uration, in reality the dust clouds are likely clumpy with some 
covering fraction, where the ionizing H II regions associated 
with the different OB asso ciations are in clos e proximity to 
the dust clouds and PDRs (IWolfire et al.lll990t) . Models fac- 
toring in the cov ering fraction prod uce a distribution of dust 
attenuation (e.g.. lDopita et al.ll2005l) . and this is more consis- 
tent with the observations. Figure|8]shows that IRAC 3.6, 4.5, 
5.8 and 8.0 fim flux densities for SMM J163554.2+661225 
are not reproduced simultaneously with the far-IR emission 
by any of the best-fitting models. Thus, spherical symmetry 
of the starburst and obscuring PDR seems insufficient to de- 
scribe both the direct stellar emission and the far-IR emission 
for this galaxy. Therefore, we conclude that star formation 
in SMM J163554. 2+661225 corresponds to star-forming re- 
gions and starbursts in local galaxies with intrinsic luminosi- 
ties of Lir ss 10 10 - 10 11 Lq, but that have been "scaled-up" 
by one-to-two orders of magnitude. Moreover, this "scaling- 
up" of star formation may b e common at high redshift. For ex- 
ample, [Ticconretly] (|200g) argued that the significantly more 
luminous sub-mm galaxies at z ~ 2 - 3 resembled "scaled- 
up" and more gas-rich versions of local ULIRGs. For SMM 
J 163554.2+66 1225, R08 found that the rest-frame mid-IR 
spectrum from 6-10 /jm is consistent with the spectra of lo- 
cal starburst galaxies and inconsistent with the spectra of local 
ULIRGs. Here, our analysis of the far-IR data, coupled with 
the SFR implied by the luminosity of the Pact line, leads us to 
a similar conclusion. 

The model we constructed in the Appendix using the 
observed rest-frame UV-near-IR SED consists of a dou- 
ble stellar population. The stellar population that domi- 
nates the bolometric emission (and the total SFR) in SMM 
J 163554.2+66 1225 is very obscured by dust that is optically 
thick to UV photons. The stellar population of the subdom- 
inant model is less obscured, and optically thin to UV pho- 
tons. This analysis is consistent with the scenario discussed 
above, where the starburst in SMM J163554.2+661225 corre- 
sponds to OB associations in close proximity to PDRs with 
a clumpy distribution. A similar configuration of multiple 
star-forming components with variable dust extinction is ob- 
served in the spatially resolved colors of other galaxies at 
z ~ 2-3, and especially those with eviden ce for recent merg- 
ers dPapovich et al.ll2005l: ILaw et al.ll2007l) . 

It remains to be seen whether the results discussed here 
for SMM J163554.2+661225 are typical of other high- 
redshift galaxies. Our full observing program will pro- 
duce similar data for eleven gravitationally lensed galax- 
ies in addition to SMM J163554.2+661225, and these will 
span a wide range of mass, optical, and IR properties. For 
SMM J163554.2+661225, the baryon (stellar + gas) mass 
we derived from the rest-frame-UV-to-near-IR SED and the 
dynam ical mass from molecular observations (IKneib et alj 
120051) both suggest M ~ 1.5 x 10 10 M . This mass is 
typical of "L*" UV-selected galaxies at these redshifts 
dPapovichet al.ll200U IShapley et aT1l2001l) . and such objects 
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likely dominate the SFR density at this redshift dReddv et alj 
120081) . However, the dust obscuration and IR luminosity are 
much larger in SMM J 163554. 2+66 1225 compared to the 
UV-selected samples, and more typical of sub-mm-sel ected 
objects (e.g. IB lain et al.1 120021: iChapman efail 120051) and 
IR-luminous ^-band-se l ected objects (e.g ., iPapovich et alj 
l2006t iDaddi et al.l [20071; IWuvts et al.1 120081) . The fact that 
the molecular- gas-emission centr oid lies between the UV- 
components led lKneib et al.l (120051) to argue that the high SFR 
in SMM J 163554.2+66 1225 results from the merger of two 
progenitors whose properties are similar to those of "typical" 
UV-selected galaxies. In this case, we are observing SMM 
J 163554. 2+66 1225 during perhaps a short-lived stage of en- 
hanced star formation. Our study of the SFR indicators in 
SMM J163554.2+661225 provides the first evidence that the 
total IR luminosity is proportional to the number of ionizing 
photons in these situations. 



5.5. Final Thoughts 

A potential source of bias is that our observations 
are sensitive to the integrated emission from SMM 
J163554. 2+661225, and we are combining datasets with a 
wide range of angular resolution. This is of concern as as 
our analysis corresponds to the flux-weighted average of the 
star-formation properties of an individual galaxy. The analy- 
sis of the rest-frame UV-to-near-IR SED (see the Appendix) 
suggests there are at least two stellar populations, with very 
different extinction properties. Nevertheless, our results prob- 
ably apply to the dominant, most luminous star-forming com- 
ponents, which dominate the nebular emission (Paa) and the 
IR emission. Therefore, we suspect that our results are valid 
for the average properties of SMM J163554.2+661225. We 
note that similar issues arise in the analysis of the integrated 
properties of the starburst galaxies (lEngelbracht et al] 120051 : 
lAlonso-Herrero et al.1 12006c ICalzetti et al.l 120071) . and some 
caution must be applied then when comparing these with, for 
example, the analysis of individual H II regions which are re- 
solved in both the narrow-band imaging of the hydrogen re- 
comb ination lines and the mid-IR data (ICalzetti et alj|2005l 
120071) . Achieving resolved observations of the Paa emis- 



sion and far-IR emission in SMM J163554.2+661225 may 
be possible with the James Webb Space Telescope and large- 
aperture sub-mm facilities. However, further work will be 
needed to quantify possible biases in the integrated emission 
from galaxies. 

6. SUMMARY 

We observed the galaxy SMM J163554.2+661225 at z = 
2.515 with Spitzer spectroscopy from 5.8-8.0 pm. These data 
are the deepest Spitzer/IRS data taken with the SL2 module of 
any galaxy of which we are aware. SMM J163554. 2+661225 
is a sub-millimeter-selected infrared-luminous galaxy, and is 
lensed gravitationally with a magnification of p = 22 by the 
rich galaxy cluster Abell 2218. This galaxy is maintaining 
a high rate of star formation based on its IR emission, and 
it has no evidence for an AGN in either its rest-frame UV 
or optical spectra, nor based on its (lack of) X-ray activity. 
We find that the rest-frame-UV-to-near-IR SED of SMM 
J 163554.2+66 1225 is best represented by the superposition 
of a double stellar population with a varying amount of dust 
attenuation, and with a total stellar mass of ~10 10 M Q (cor- 
rected for the gravitational lensing magnification). In this 



model the stellar population that dominates the star formation 
is heavily extincted, A(V) ~ 3 mag, while the subdominant 
stellar population is optically thin to UV photons. 

We detected the Paa emission line in the Spitzer spectrum 
with a redshift z = 2.515 ± 0.003. The luminosity of the Paa 
emission line is L(Paa) = (2.05 ±0.33) x 10 42 erg s" 1 , cor- 
rected for the gravitational-lensing magnification, with a rest- 
frame equivalent width EWq = 363 ± 56 A. We compared the 
Paa luminosity to the Ha luminosity and derived a nebular 
extinction of A(V) = 3.6 ± 0.4 mag, although this depends on 
the assumed attenuation law and remains a systematic un- 
certainty. This is consistent with that derived from model- 
ing of the galaxy's rest-frame UV-to-near-IR SED, and im- 
plies that the dust attenuation affects the ionized gas and stars 
uniformly. The dust-corrected Paa luminosity is L(Paa) = 
(2.57 ±0.43) x 10 42 erg s" 1 , corrected for the gravitational 
lensing magnification. This corresponds to an ionization rate 
of go = 1-6 x 10 55 7 s~\ implying SMM J163554.2+661225 
contains on the order of 10 6 O stars. Assuming an IMF 
with a Salpeter-like slope from 0.1 to 100 M Q yields a SFR 
ip= 171±28M yr" 1 . 

The Paa-derived SFR agrees with the upper range of SFR 
implied by the total IR luminosity. The uncertainty on 
the total IR luminosity is a factor of order two, primarily 
due to the lack of data at ^100-300 pm (rest-frame ~30- 
100 pm). The measured Paa and rest-frame 8 fim luminosi- 
ties are consistent with the extrapolated relation observed in 
local galaxies and star-forming regions. This implies that 
both the monochromatic 8 pm and total IR luminosities are 
dominated by heating from ongoing star-formation and they 
are proportional to the number of ionizing photons. How- 
ever, the measured rest-frame 24 pm luminosity is signif- 
icantly lower in SMM J 163554.2+66 1225 compared to lo- 
cal galaxies with comparable Paa luminosity. Thus, SMM 
J 163554.2+66 1225 appears to lack a warmer dust component 
{To ~ 70 K), which is typical in local galaxies of comparable 
L(Paa). The nature of SMM J163554.2+661225 seems very 
different from the properties of local ULIRGs, even though 
they have comparable SFRs. Comparing the IR emission of 
SMM J163554. 2+661225 to expectations from empirical and 
radiative transfer models, we conclude that the starburst in 
SMM J 163554.2+66 1225 has similar physics as those in lo- 
cal galaxies with intrinsic luminosities, ~ 10 10 — 10 11 L@, but 
"scaled-up" by one-to-two orders of magnitude. 

The implied stellar and dynamical masses of SMM 
J163554.2+661225 are consistent with those of typical "L*" 
UV-selected objects, which dominate the SFR density at this 
redshift. The implied timescales and starburst age are on the 
order of the dynamical time, implying that this galaxy is only 
part way through its elevated star-formation episode. Our 
analysis here provides the best measurement yet of the SFR 
in a galaxy involved in such an episode and it shows that the 
IR-luminosity traces the total SFR in this situation. While 
our analysis here pertains to only one galaxy of any signif- 
icant redshift, as we extend this work to our larger sample, 
it will allow us to study galaxies spanning a wider range of 
star-forming properties and luminosities. 
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APPENDIX 
STELLAR POPULATION FITTING 

Photometry is available for SMM J 163554.2+66 1225 from 
HST/WFPC2 (BVI), WHT/INGRID (JK S ), and Spitzer/IRAC 
(channels 1-4), spanning 0.5-8 /jm in the observed frame 
(K04, R08). We show the SED of this galaxy in figure|9] We 
fitted a suite of stellar population synthesis models to these 
data to estimate the properties of th e stellar populations in this 
galaxy. We used both the models o flBruzual & Chariot! (120031) 
and the updated 2007 version, which include a more physi- 
cally motivated treatment of the thermally pulsating asymp- 
totic giant branch stars. We use models with a Chabrier IMF, 
although for consistency with other SFR indicators, we mul- 
tiply the derived stellar masses and SFRs by a factor of 1.8 
to convert them to equivalent quantities using a Salpeter-like 
IMF. We allow for a range of exponentially declining star- 
formation histories with e-folding time r = 1 Myr ("instan- 
taneous burst") to r = 100 Gyr ("constant" star-formation), 
metallicities of 0.2 and 1.0 Z w , and we a llow a range of dust 
extinction using the lCalzetti et alJ (120001) law with E(B-V) = 
- 1 . We note that this dust law may be inappropriate for 
hea vily obscured sources such as SMM J 163554.2+66 1225 
(see iGoldader et ani20()2). but it provides a useful compari- 
son to other work modeling galaxy spec tral energy distribu- 
tions (SEDs). See IPapovich et all (hooil 120061) for details of 
the SED fitting. 

We find that the photometry are not well fit by simple 
models with single component, exponentially declining star- 
formation histories, as illustrated in figure [9] Moreover, the 
best-fit single-component model is unphysical. It has a stellar 
mass of, M = 7 x 10 10 M Q (corrected for the gravitational lens- 
ing magnification), stellar population age, t = 19 Gyr (greatly 
exceeding the age of the Universe at this redshift) , formed 
with a star-formation e-folding timescale, r = 200 Myr, and 
with a color excess, E(B-V) = 0.32. The reduced % 2 for the 
best-fitting single- component model is x 2 / v = 14- 

SMM J 163554.2+66 1225 likely consists of multiple star- 
forming components, with variable extinction and star- 
formation properties, all of which contribute to the rest-frame 
UV through IR emissio n. Such situations are seen in local IR- 
lumin ous galaxies (e.g.. lCharmandaris. Le Floc'h. & Mirabell 
120041) . and arguably apply to high redshift galaxies as well. 
As discussed by K04, HST and ground-based imaging shows 
that SMM J163554.2+661225 consists of several distinct blue 
components, with a luminous red core. We therefore tested 




Observed Wavelength [yu.m] 

FIG. 9.— Spectral energy distribution of SMM J163554.2+661225 from 
0.5 to 8 fim in the observed frame. The data points show the measured 
photometry taken from K04 and R08, uncorrected for gravitational lensing. 
We find that simple, single-component models with exponentially declining 
star-formation histories do not reproduce the photometry well. The dashed 
line shows the best-fit single-component model, with parameters given in 
the text. However, a double stellar-population model with two star-forming 
components better reproduces the photometry, which is illustrated by the 
solid curves in the figure. The red curve corresponds to the "very obscured 
component" and the blue curve to the "less obscured component", with model 
parameters given in the text. The thick black curve shows the sum of these 
two components. 

SED models consisting of two star-forming components. We 
find that this "Double Stellar Population" produces a better 
fit to the photometry, where the two components correspond 
to a very dust-obscured star-forming component, and a less 
obscured component. Figure [9] shows a characteristic model 
that reproduce the data. The reduced \ 2 of the fit to the 
data for the models in the figure is \ 2 l v — 4-0, greatly im- 
proved over the single-component model above. In these 
models, the "very obscured component" dominates both the 
SFR and the total stellar mass. In the figure, the very ob- 
scured component consists of a stellar population forming 
with a constant SFR with an age, t = 80 Myr, and a color 
excess, E{B-V) = 0.80, corresponding to A(V) = 3.2 mag, 
stellar mass, M= 1 x 10 10 M , and SFR, ip = 110 M Q yr" 1 
(the stellar mass and SFR have been corrected for the grav- 
itational lensing magnification). The derived color excess is 
much larger than those derived fo r UV-selected galax ies at 
this redshift (IPapovich et alj|2001b IShaplev et alj(2001l) , but 
is approximately in the range of those der ived for sub-mm 
selected at this redshift dBorvs et alJ 120051) . Also, the dust 
extinction we derived from modeling the rest-frame UV-to- 
near-IR SED here is consistent with that derived from the Ha 
and Paa emission-line ratio, and implies that the dust attenua- 
tion affects the ionized gas and stars uniformly. The subdom- 
inant "less obscured component" of the double-component 
model consists of a young stellar population, t = 40 Myr, with 
a star-formation e-folding timescale, t = 50 Myr, moderate 
color excess, E(B-V) = 0.30, stellar mass, M = 1 x 10 8 M©, 
and SFR, ip = 3 M© yr" 1 (the stellar mass and SFR have 
been corrected for the gravitational lensing magnification). 
The extinction in this component is c onsistent with the up- 
per range seen in L BGs at this redshift ( Pap ovich et alj|200lt 
IShaplev et alJl200lh . 

Combined with the total gas mass, M gas = 4.5 x 10 9 M© 
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dKneib et all 120051) . the total baryonic mass is M* w 1.5 x 
10 10 Mq, which is consistent with the estimate of the dynam- 
ical mass of the sy stem as inferred from the molecular gas 
dKneib et al.l 120051) . Therefore, we conclude generally that 
while the SED of SMM J 163554.2+66 1225 is complex, it is 



broadly consistent with a multiple component starburst, with 
total stellar mass ~ 10 10 M Q , where a deeply obscured 
component is responsible for for the vast majority ( >80%) of 
the SFR. 
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